Photoionization in negative streamers: fast computations and two propagation modes 
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Streamer discharges play a central role in electric breakdown of matter in pulsed electric fields, 
both in nature and in technology. Reliable and fast computations of the minimal model for negative 
streamers in simple gases like nitrogen have recently been developed. However, photoionization was 
not included; it is important in air and poses a major numerical challenge. We here introduce a fast 
and reliable method to include photoionization into our numerical scheme with adaptice grids, and 
we discuss its importance for negative streamers. In particular, we identify different propagation 
regimes where photoionization does or does not play a role. 

PACS numbers: 52.80.Mg, 52.27. Aj, 52.65.Kj 



Streamers are a generic initial stage of sparks, light- 
ning and various other technical or natural discharges . 
More precisely, when a high voltage pulse is applied to a 
gap of insulating matter, conducting streamer channels 
grow through the gap. Streamer propagation is charac- 
terized by a strong field enhancement at the channel tip. 
This field enhancement is created by a thin curved space 
charge layer around the streamer tip as many computa- 
tions show. Such computations are quite challenging due 
to the multiple inherent scales of the process. 

Recent streamer research largely concentrates on pos- 
itive streamers in air or other complex gases for indus- 
trial applications This is because positive streamers 
emerge from needle or wire electrodes at lower voltages 
than negative ones . Natural discharges such as sprites 
0, on the other hand, occur in both polarities 0, in 
particular, when they are not attached to an electrode 
and therefore double ended. Photoionization (or alter- 
natively background ionization) is essential for positive 
streamers: as their tips propagate several orders of mag- 
nitude faster than positive ions drift in the local field, a 
nonlocal photon-mediated ionization reaction is thought 
to cause the fast propagation of the positive ionization 
front. Negative streamers, on the other hand, have ve- 
locities comparable to the drift velocity of electrons in 
the local field, therefore a local impact ionization reac- 
tion can be sufficient to explain their propagation. This 
is why photoionization in negative streamers has received 
much less attention, most recent work concentrating on 
sprite conditions with relatively low electric fields |5| . 

The nonlocal photoionization reaction depends 
strongly on gas composition and pressure 0, in partic- 
ular, it is much more efficient in air than in pure gases. 
Furthermore, in air its relative importance saturates 
for pressures well below 60 Torr (« 0.1 bar), while it is 
suppressed like « 60 Torr/p at atmospheric pressure and 
above. In this paper we study the effects of photoioniza- 
tion on the propagation of negative streamers by means 
of efficient computations with adaptive grids. 

Streamer model. Streamer models always contain 



electron drift and diffusion, space charge effects and the 
generation of electron ion pairs by essentially local impact 
ionization. We will use a fluid model in local field approx- 
imation as described, e.g., in Refs. A numerical 
code with adaptive grid refinement was introduced in 
to investigate negative streamers in pure nitrogen, where 
photoionization plays a negligible role. With this code 
even streamer branching could be determined accurately 
1^]. On the other hand, in gases like air where photoion- 
ization cannot be neglected, photons emitted from ex- 
cited molecules can act as a non-local source of electron- 
ion pairs; this has to be included in the computations. 
The challenge lies in maintaining computational speed 
and accuracy while introducing the nonlocal interaction. 

More precisely, the number of photoionization events 
at a given point r results from integrating the emission 
/(r') of photons at every point r' of the gas volume mul- 
tiplied by a kernel that contains an absorption function 
and a geometrical factor. The production of photons in 
air is, on the other hand, proportional to the number of 
impacts of free electrons on nitrogen molecules and hence 
can be related to the impact ionization Si{r). Thus, 



P+Pq 



(1) 



where ^ is a proportionality factor that weakly depends 
on the local reduced electric field although it is com- 
monly assumed to be constant and about ^ — 0.02. We 
must note here that, since the only data accessible from 
macroscopic observations is the product ^f{r), this is of- 
ten packed into a single function and called, by a slight 
abuse of terminology, absorption function. In this letter, 
however, we prefer to apply this term only to f{r). The 
factor pg/ {p + Pq) accounts for the probability of quench- 
ing, i.e. for the non-radiative deexcitation of a nitrogen 
molecule due to the collision with another molecule. The 
pressure Pq is called quenching pressure and will be taken 
here as pq = 60 Torr Q . There is some uncertainty over 
this value and some authors (sl lid, prefer pq — 30 
Torr. However, different values of pq within this range 
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FIG. 1: The function ^/ of the photoionizing radiation in 
the range 980 - 1025 A in air, taken from Ref. (squares) 
and fit according to l|21l (sohd hne) with the parameters Ai — 
6.0 • lO"'^ cm-^Torr~\ A2 = 3.55 • 10"'' cm"^Torr"\ Ai = 
0.059 cm"^Torr"\ A2 = 0.010 cm"^Torr~\ 

affect our quantitative results only marginally and our 
numerical approach and qualitative observations remain 
unchanged. 

Evaluating the integral numerically in each time 
step is very time consuming, since for each grid point 
r one has to add the contributions of all emitting grid 
points r'. Kulikovsky has assumed cylindrical sym- 
metry and has considered only a relatively small number 
of uniformly emitting rings, interpolating at finer levels. 
This approximation ignores the small-scale details of the 
density and electric field distributions that matter, e.g., 
in a branching event. 

Numerical implementation of photoionization. 
We here present a different numerical method that allows 
us to keep calculating with a locally appropriately refined 
numerical grid, and nevertheless to obtain reliable results 
within decent computing times. Our approach relies on 
approximating the absorption function as 

/(|r-r'|) = ^-^5]V-'^l— 'I (2) 

^ i=i 

where Ai . . . Aat and Ai . . . Aj^ fit the experimental data 
as closely as possible. This form has the advantage that 
the integral can be expressed by a set of Helmholtz 
differential equations for the Sphj as 

N 

Sph = — ~ ^ ^ ^j^phjj ^ ^"j) ^ph,j = Si, (3) 

P + Pl , = 1 

with the boundary condition Sphjiv) — > far away from 
the high field areas. Thus one now can use the very fast 
algorithms available for solving elliptic partial differential 
equations with separable variables, such as described in 
Ref. and implemented in the freely downloadable li- 
brary FISHPACK. The same algorithm was used in Ref.!3 
to solve the electrostatic problem [l^ . 

For nitrogen-oxygen mixtures like air, the most reliable 
model for / is provided by Ref. Abased on the experimen- 
tal measures of Ref. Q, despite some recent controversy 



over these data [lolll5| . Fig.^shows the data for / from 
Ref. 13 together with our fit of form (0) with N = 2. 

Note that the asymptotic behavior of (O for |r— r'| — > 
and |r — r'l 00 disagrees with that predicted by Ref.0. 
Nevertheless, these differences cannot be seen in Fig. ^ 
For very small distances between the emitting excited 
state and the ionized molecule, the impact ionization is 
dominant anyway. At distances much larger than the 
largest absorption length l/(Ajp), where most radiation 
is absorbed, the identical exponential decay in r domi- 
nates over the different powers of r. 

Similarity laws. Without photoionization, there are 
similarity laws between streamers at different pressures: 
They are equal after rescaling lengths, times and fields 
with appropriate powers of the pressure — this gener- 
alizes Townsend's historical finding that the ratio of elec- 
tric field over pressure E/p is the physically determining 
quantity in a discharge, not E and p separately. Pho- 
toionization introduces a nontrivial pressure dependence 
through the factor pq/ (p + pg) in |^ and thus breaks the 
similarity laws between streamers at ground level and 
those in the high altitude, low pressure regions where 
sprites appear l5l|. 

simulation setup. We have incorporated photoion- 
ization into the numerical code of Ref. as described 
above. Air was approximated as an oxygene-nitrogen 
mixture in the ratio 20:80. In order to study the 
effect of photoionization on streamer propagation at 
different background electric fields, we used fields of 
100 {p/po) kV/cm and of 40 {p/po) kV/cm where po is at- 
mospheric pressure. Furthermore, we studied three pres- 
sure regimes, namely atmospheric pressure (760 Torr) 
and 0.05 Torr, which corresponds to the pressure of 
the atmosphere at around 70 km above sea level, where 
sprites are commonly observed, and also the case without 
any photoionization, which corresponds to the physical 
limit of very high pressures, when all excited states are 
rapidly quenched, or to the case of pure nitrogen. 

The length of the computational domain was 
4.7mm/(p/po) for the higher and 9.4mm/(p/po) for the 
lower electric field. The radial extension was large enough 
that the lateral boundaries did not influence phenom- 
ena. At the cathode we imposed homogeneous Neu- 
mann boundary conditions, roughly equivalent to a free 
electron inflow into the system. An initial ionization 
seed was introduced near the cathode as an identical 
Gaussian density distribution for electrons and ions with 
a maximum of 8.2 • 10^ / (p/po)^ imn~^ and a radius of 
23(po/p) /im. 

Simulation results and conclusions. Some simu- 
lation results for the evolution of the streamer head in 
different fields and pressures are shown in Fig. 13 Let us 
focus first on the high field regime which is represented 
in the left column of the figure; there it can be seen that 
during the first three to four time steps, the streamer 
development is barely affected by photoionization pro- 
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FIG. 2: Evolution of streamers in a field of f = 100 kV/cm 
p/po (left column) and in a field of f = 40 kV/cm p/po (other 
three columns). Plotted is the contour of the half-maximum 
of the space charge at different times. The time interval be- 
tween two consecutive snapshots is At = 150 ps / (p/po) for 
the high field and At = 2400 ps / (p/po) for the low field. 
Lengths are measured in units of 2.3 / (p/po). Shown 
are streamers (a) without photoionization, (b) in air at at- 
mospheric pressure, and (c) at low pressures (p <g 60Torr). 
Note that the computational domain is larger than the plotted 



cesses. However, eventually a new phase sets in where 
the streamer accelerates significantly. This acceleration 
is the stronger, the higher the relative contribution of 
photoionization, i.e., the lower the pressure. On the other 
hand, field enhancement is much weaker: it increases by 
~ 400% without photoionization and only by ^ 60% in 
the low pressure case. 

In the lower field case, a very different behavior is 
seen: photoionization hardly changes the streamer veloc- 
ity. However, now it does suppress streamer branching 
as also found in Ref. 0. This can directly be related to 
the fact that photoionization makes particle distributions 
smoother, and that a smoother space charge layer is less 
susceptable to a Laplacian instability 

This smoothening dynamics can be made more precise 
by plotting the logarithm of the electron density along 
the symmetry axis of the streamer in Fig. 13 Photoion- 




FIG. 3: Evolution of the logarithmic electron density on the 
streamer axis at pressure p <C 60 Torr (solid lines) and with- 
out photoionization (dashed lines) in the high field (above), 
corresponding to columns 2 and 4 in Fig. |5| and low field (be- 
low), which corresponds to the upper and lower plots in the 
left column of Fig. |21 



ization creates a smoothly decaying density tail ahead of 
the ionization front that initially is not visible on a lin- 
ear (non-logarithmic) scale. The point where the steep 
density decrease crosses over a smoother photoionization 
induced decay, moves toward higher density levels with 
time. For low fields (Fig. O below) up to the time when 
the streamer without photoionization branches, the large 
density levels visible in Fig. |21 move essentially with the 
same velocity. This is different in the high field case 
(Fig. 13 above): there the photon created leading edge 
eventually dominates the complete decay of the electron 
density and pulls the ionization front to much higher ve- 
locities 17J. 
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